Generation of High Brightness Electron Beams via Ionization Induced Injection by 
Transverse Colliding Lasers in a Beam-Driven Plasma Wakefield Accelerator 
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The production of ultra-bright electron bunches using ionization injection triggered by two trans- 
versely colliding laser pulses inside a beam-driven plasma wake is examined via three-dimensional 
(3D) particle-in-cell (PIC) simulations. The relatively low intensity lasers are polarized along the 
wake axis and overlap with the wake for a very short time. The result is that the residual momentum 
of the ionized electrons in the transverse plane of the wake is much reduced and the injection is 
localized along the propagation axis of the wake. This minimizes both the initial thermal emittance 
and the emittance growth due to transverse phase mixing. 3D PIC simulations show that ultra-short 
(~8 fs) high-current (0.4 kA) electron bunches with a normalized emittance of 8.5 and 6 nm in the 



two planes respectively and a brightness greater than 1.7 x 10 19 A • rad 2 
realistic parameters. 



m can be obtained for 



PACS numbers: 52.38.Kd, 41.75.Jv, 52.35.Mw 

The demonstration of the Linac Coherent Light Source 
(LCLS) as an X-ray free electron laser (X-FEL) [l| has 
given impetus to research on the fifth-generation light 
sources [2|. The goal is to make X-FELs smaller and 
cheaper while decreasing their wavelength and increas- 
ing their coherence and intensity. The FEL performance 
is partially determined by the brightness of the electron 
beam that traverses the undulator. The brightness is 
defined as B n = 2//e 2 where / is the beam current 
and e n is the normalized emittance of the beam. In or- 
der to make the length of the undulator needed to drive 
the SASE-FEL ^ into saturation, shorter, high current 
(^kA), multi GeV electron beams with e n ~ lOnra will 
be needed. These emittances are an order of magnitude 
smaller than those from state-of-the-art photoinjector RF 
guns [4]. In this letter, we show the generation of ultra- 
bright electron bunches using ionization injection trig- 
gered by two transversely overlapping laser pulses inside 
a beam-driven wake in plasma. In our scheme, the rel- 
atively low intensity lasers are polarized along the wake 
axis and overlap with the wake for a very short time. 
Particle-in-cell (PIC) simulations using OSIRIS ^ show 
that this geometry reduces the residual momentum of 
the ionized electrons in the transverse plane and local- 
izes them along the propagation axis of the wake leading 
to an electron beam with a brightness greater than 10 19 
A-rad -2 -m -2 that would be highly attractive to future 
light sources. 

When a dense (n^ > n p , k v a r ^ z < 1), ultra-relativistic 
(7 ^> 1) electron beam propagates through a plasma, 
the plasma electrons can be completely blown out by the 
repulsive Coulomb force of the beam leaving behind a 
cavity of more massive ions [6-8] which then pull the 
electrons back creating a wakefield with a phase velocity 



equal to the velocity of the beam. Here n^, n p , k p and <j rjZ 
are beam density, plasma density, inverse of the plasma 
skin-depth and transverse and longitudinal r.m.s. size of 
the electron beam respectively. The accelerating and fo- 
cusing fields inside this wakefield have ideal properties for 
acceleration of electrons while maintaining beam quality 
[(Hsl and high-gradient acceleration by such wakes has 
been experimentally demonstrated [9 12]. 

For a plasma density ~ 10 18 cm -3 , the wavelength 
of the ion cavity is about several tens of microns mak- 
ing synchronization and efficient capture of externally in- 
jected electrons into such a cavity extremely challenging. 
Self-injection of electrons in plasma wakes is conceptu- 
ally simple, however, it still can not generate sufficiently 
high brightness beams needed for next generation light 
sources [lH [3 • Other electron injection schemes such as 
ponderomotive force injection have been proposed [l5[ 
for laser-plasma wakefield accelerators. Here a pump 
laser pulse excites the wake and a second injector pulse 
propagates orthogonally to the pump disrupting the or- 
bits of some of the plasma electrons leading to trapping. 
Soon thereafter the collinear colliding pulse injection was 
suggested [l6j and experimentally demonstrated 17]. In 
addition a sudden [l8j or gradual [l9| density transition 
from a high plasma density to a low plasma density has 
also been shown to inject particles into plasma wakes. 
Another technique is ionization injection where electrons 
are produced inside the wake by the electric field of a 
laser pulse or the drive electron beam (that produce the 
wake) where they can be more easily captured and accel- 
erated. Ionization injection is attractive because it offers 
the potential to control the charge and emittance of the 
accelerated beam. Very recently it was proposed to com- 
bine the ionization injection via an auxiliary laser pulse 
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into a beam driven wake [20j. This approach allows the 
use of a lower intensity ionizing laser, thereby further 
reducing residual momentum and hence the transverse 
emittance of the injected electrons. In this letter, we 
show that electron injection into a beam-driven PWFA 
via tunnel ionization in the overlap region of two laser 
pulses (moving transversely across the wake) can gen- 
erate an electron beam with extremely small transverse 
emittances and therefore an extremely high brightness. 

This mechanism is explored using the 3D PIC code 
OSIRIS [5[ in Cartesian coordinates using a moving win- 
dow. We define the z-axis to be the propagating direc- 
tion of the drive beam, and the x-axis to be the propa- 
gating direction of the colliding laser pulses with their 
electric field polarized along the z-axis. The simula- 
tion window had a dimension of 89 x 81 x 121 /im with 
1400 x 512 x 760 cells in the x, y and z directions respec- 
tively. This corresponds to cell sizes of 0.5/c^" 1 in the x 
direction and 1.25/c^" 1 in the y and z directions. The code 
uses Ammosov-Delone-Krainov (ADK) tunneling ioniza- 
tion model [211 ]. 

For simplicity, the simulation is initialized with plasma 
with a density of n p = 2.4 x 10 17 cm -3 , represented 
by 8 particles per cell, and neutral He with a density 
of 1.1 x 10 18 cm -3 represented by 8 neutral atoms per 
cell. The pre-ionized plasma can be viewed as a fully 
ionized separate gas. A 500 MeV drive beam of the form 
_^_ ?e - 2 /2^ e -z 2 /2^ with Gz = 11.4 ^ m> ar = 7 . 6 

/im, and N = 1.25 x 10 9 (200pC) respectively, propa- 
gates through the plasma and excites the wake. The self 
electric field of the beam (^50 GeV/m) does not ionize 
the helium atoms. In addition, two counter propagating 
laser pulses moving along the + and - x-axis directions 
are synchronized with the electron beam so that they 
overlap inside the first bucket (ion cavity) near the point 
where the longitudinal electric field vanishes, z.e., where 
E z = 0. Each laser has a normalized vector potential of 
do = 0.016, a duration of r = 20 fs, and a focal spot 
size of wo = 6 /im. These parameters correspond to each 
laser having a focused intensity of 5.5 x 10 14 W/cm 2 for 
a wavelength of 800 nm. 

We first examine the injection process. It is easier 
to trap and to control the self-injection of an electron 
that is born (ionized) at rest inside the wake. The trap- 
ping threshold can be written as [22| Aip = ip — -0 init < 
— 1 + (1 + [pi/mc] 2 ) 1 ^ 2 h P h where if) = e{(j> — A z )/mc 2 is 
the normalized wake potential and ^i n it is the wake po- 
tential when and where an electron is created. The wake 
potential is at a minimum at the rear of the cavity so it is 
easiest to trap an electron if it is created when the poten- 
tial is at a maximum (a zero for E z ) which occurs in the 
middle of the cavity. Furthermore, low emittance beams 
are generated if the electrons are born near the trans- 
verse axis. For laser ionization, an electron is born inside 
a laser so it also acquires a residual drift in the polariza- 
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FIG. 1. Snapshots from PIC code simulations illustrating the 
transverse colliding pulse injection of helium electrons into 
the ion cavity. Snapshots (a) to (c) show the charge density 
distribution of driver beam, wake electrons and helium elec- 
trons at three different times (a) ^80 fs before laser pulses 
collision (b) around laser pulses' collision time (c) ~200 fs 
after collision when the injected electrons become trapped in 
the wake. The ionized electrons seen on the left of each figure 
are due to the wakefields themselves exceeding the ionization 
threshold for helium at the rear of the first bubble however 
none of these electrons are trapped and accelerated in the 
wakefield. 

tion direction of the laser of p/mc = eA in \ t /mc 2 (approx- 
imately in multi-dimensions). The amount of charge that 
will be injected is controlled by the neutral gas density. 

Fig. [1] illustrates the injection process through a se- 
ries of snapshots that shows the formation of the wake, 
the collision of the two lasers triggering injection within 
this wake through ionization, and the trapping of the in- 
jected electrons, where the densities of the drive electron 
beam, wakefield electrons and laser-ionized electrons are 
plotted with colors black, blue and pink, respectively. In 
each snapshot a stable highly nonlinear wakefield cavity 
is seen. In Fig. [TJa) the injector laser pulses are seen just 
moving across the sharp electron sheath of the ion cav- 
ity. Because of the low intensity of laser pulses, approx- 
imately 5.5 x 10 14 W/cm 2 , the ponderomotive force of 
the laser pulse is not enough to perturb the sheath elec- 
trons, and thus the wake structure remains unaffected. 
In Fig. DIb), the lasers collide on the axis, where they 
have the maximum (overlapping) intensity at the posi- 
tion where E z = 0. The local laser intensity exceeds the 
ionization threshold only where the lasers overlap (each 
laser has an intensity below the ionization threshold) and 
a large fraction of neutral helium atoms within this vol- 
ume is now ionized. As the laser pulses travel past the 
collision point the injection ceases. Note that the drift 
momentum of the electrons is along the laser polariza- 
tion direction and therefore now affects predominantly 
the longitudinal momentum spread of the beam, leading 
to a longitudinal emittance of 0.06 keV-ps in this exam- 
ple. These laser ionized helium electrons then respond to 
the wake fields and are rapidly accelerated to a longitu- 
dinal velocity close to the phase velocity of wakefield by 
the time they have slipped backwards to the rear of the 
ion cavity. They then begin to move synchronously with 
the wake, as depicted in Fig. [He). 
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FIG. 2. Ion density (color-coded) indicates the ionization 
level of helium gas. Wake electrons are colored blue, (a) 
the ionization level at the instant of collision, (b) 80 fs after 
the collision, when the laser pulses are on the verge of exiting 
the electron sheath, (c) The lineout on the red dashed line 
in (a). Laser- ionized He + is shown in blue, sheath electron 
density is shown in green and black dotted line indicates the 
initial ion background. 



The details of the ionization process can be more 
clearly seen by plotting the density of the He + ions. The 
superposition of the two lasers gives rise to a standing 
wave with a wavelength Ao/2, a frequency cjo, and an 
amplitude four times that of a single pulse. Fig. [2] shows 
the He + ion density distribution resulting from laser ion- 
ization at two different times. The laser created He + 
ions are immobile because they are more massive, thus 
their location represents the birthplace of the electrons. 
In Fig. Efa), one can observe the layered ion distribu- 
tion reflecting the standing wave of electric field of the 
two laser pulses at the instant where the lasers overlap. 
In Fig. [2fc) a lineout corresponding to the red dashed 
line of Fig. [2ja) is shown. The inset shows that under 
such laser intensity, He atoms are almost fully ionized to 
He + at the antinodes of the overlapping intensity. He- 
lium electrons are mainly born within 2 /im of the axis. 
In Fig. [2fb) we show a snapshot of ion density 80 fs af- 
ter the laser collision at which time the laser is near the 
sheath. The superposition of the laser and wake fields 
near the electron sheath leads to some off-axis ioniza- 
tion. Additionally, there is ionization from the wakefield 
alone near the rear of the bubble. Fortunately, in this 
case these off-axis electrons (and those in the rear) do 
not get trapped because their initial positions are too 
close to the rear and too far from the axis of the cavity. 
However, for other beam and plasma parameters, these 
undesired He electrons can get trapped deteriorating the 
beam quality. 

Fig. 3 illustrates the phase space distributions in 
each plane about 500 fs after the collision. The 




FIG. 3. The (a) x-p x , (b) y-p y and (c) z-p z phase space 
distribution about 500 fs after the pulse collision. 



projected normalized emittance is obtained by e n = 
Pi V '{%' 2 ) (% 2 } — {xx') 2 for each of the two transverse 
planes. Initially, the injected beam has an ultra- low 
projected transverse normalized emittance for the whole 
bunch, of about 8.5 nm in the x direction and 6 nm in 
the y direction, which is observed to be invariant after 
propagating 100 /im. At this distance the beam has an 
average energy 5.3 MeV, a slice energy spread of ~12 keV 
and a total charge of 4.6 pC. The projected longitudinal 
emittance as explained earlier is also very small, about 
0.06 keV-ps. The beam current profile is near flattop, 
with a r.m.s. pulse duration around 8 fs and peak cur- 
rent 0.44 kA. The brightness of the beam B n — 21 v je 1 n is 
estimated to be 1.7 x 10 19 A • rad -2 • m -2 , more than 3 or- 
ders of magnitude higher than that of LCLS. In principle, 
the trapped beam charge and the normalized emittance 
can be further optimized by adjusting the neutral gas 
concentration and the spatial sizes of injector pulses. 

To achieve the extremely small e n , certain condi- 
tions have to be fulfilled. First, the electrons have to 
"start" with small transverse momenta. In the proposed 
scheme, since the laser pulses are polarized along z axis 
and propagate perpendicularly to the x axis, the elec- 
tric (and vector) field components E Xl E yi E z scale as 
E x ~ € 2 E Z , E y ~ eE z , where e is a characteristic small 
parameter defined as e = 1/fco^o? around 0.02 in our 
simulation, therefore the residual transverse momenta 
Px->Py of electrons just after ionization scale as p x /m e c ~ 
e 2 ao, p y /m e c ~ eao, which are extremely small. Second, 
the transverse electron beam sizes are determined 
by the laser intensity contour above ionization thresh- 
old near collision point, approximately oc cr, which 
for the simulated parameters is around a few microns. 
Combining the above two factors, the intrinsic or "ther- 
mal" emittance of the electron beam, which is defined as 
eth = (p±)wb, is very small. 

However, in order to preserve the "thermal" emittance, 
it is necessary to avoid phase mixing , which arises when 
electrons within the bunch are born at different times and 
therefore at different phases of their betatron oscillations 
[2^ . To illustrate this point, using OSIRIS we simulate 
the injection process of both the transverse colliding pulse 
injection scheme (this paper) and the recently proposed 
longitudinal injection scheme [lo[. In the longitudinal in- 
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FIG. 4. Comparison of x-p x phase space evolution between 
the transverse colliding pulses injection (left column) and the 
longitudinal injection (right column) at the instant 130 fs 
(top) and 260 fs (bottom) after the onset of injection. Phase- 
mixing occurs during the longitudinal injection. 



jection scheme, a single injection laser pulse propagates 
colinearly at an optimum distance (67 fs) behind the 
beam driver, also at the phase where E z = 0. For both 
simulations, the drive beam is the same, with 200 pC of 
charge and a r = a z = 3.8 /im, and the pre- ionized plasma 
density is set to 5 x 10 1T cm -3 . The simulation domain is 
63.5/im x 50.8/im x 63.5/im and 1000 x 400 x 500 cells were 
used in the transverse injection case and 500 x 400 x 1000 
cells in the longitudinal injection case. The injection 
lasers have the same pulse duration (r = 20 fs) and a 
FWHM spot size (5 /im), but with different focal inten- 
sities. For longitudinal injection, the laser is polarized in 
x direction and focused with ao = 0.035 or an intensity 
of 2.6 x 10 15 W/cm 2 . For the colliding pulse injection, 
each laser has a ao = 0.016, therefore near the antin- 
ode of the standing wave, the combined laser intensity is 
just above the ionization threshold of helium. The neu- 
tral helium densities are also different. The He density 
is 1.1 x 10 18 cm -3 , and 5.2 x 10 16 cm" 3 for the "trans- 
verse" and " longitudinal" propagation cases respectively. 
By doing so, both cases obtain similar injected charge 
(1.9 pC for transverse injection and 3 pC for longitudi- 
nal injection). If we choose much higher helium density 
in the longitudinal injection case, the injected charge can 
be larger but the emittance gets much worse, e.g., for the 
same He density as that used in the transverse injection 
case e n was ten times worse. We also have simulated the 
colliding pulse geometry with the lasers polarized along 
y instead of x. In this case the thermal emittance along 
y is comparable to that of the longitudinal case illus- 
trating the importance of minimizing both the thermal 
emittance and phase mixing. 

Fig. H] illustrates that the x-p x phase space evolution is 
very different between the transverse colliding pulses in- 



jection (left column) and the longitudinal injection (right 
column). Fig. |4{a) and (b) are taken 130 fs after the on- 
set of injection (at this time injection has ceased for col- 
liding pulses while it is still ongoing for the longitudinal 
injection scheme). Fig. H^c) and (d) are taken about 260 
fs after the onset of injection. There is phase space rota- 
tion in the colliding pulse case (as expected) and there is 
far less phase space mixing as compared to longitudinal 
injection, because the injection distance in the longitudi- 
nal scheme is much longer than that in the colliding pulse 
scheme. The injection distance in the longitudinal case 
is on the same order as the Rayleigh length of the laser 
pulse, zr = ttWq/Xo. Over this distance He electrons 
born at different times have different Betatron motion 
phase <£. We define A<£> as the betatron phase differ- 
ence between the first and the last ionized electron. As 
seen in Fig. H^b), the first ionized electrons have rotated 
to the blue dashed line, while the final ionized electrons 
have just been released (roughly in the red dashed box). 
Xu. et al. have studied the electron dynamics of ioniza- 
tion injection and shown that the phase divergence A<£ is 
mainly determined by injection distance A as A<1> oc \/~K 
[23j| . The implication being that a longer injection dis- 
tance will cause more phase mixing leading to a larger 
final emittance. Therefore, the effect of phase-mixing is 
found to be small in the proposed scheme, because all the 
He electrons are released in short time which is controlled 
by the overlap time of laser pulses. As shown in Fig. SJ^a) 
and (c), the electrons rotate in the phase space with a 
small phase divergence, similar to a laminar beam. In 
our simulations, the normalized emittances in the longi- 
tudinal injection scheme are 34 nm in x direction and 25 
nm in y direction, both four times larger than those ob- 
tained in the transverse colliding pulses injection scheme. 
We note that if an energy chirp grows along the beam, 
then different slices of the beam will oscillate at different 
betatron frequencies leading to a growth in the projected 
(not slice) emittance. 

We have examined the effect of power imbalance be- 
tween the two colliding laser pulses on the injection pro- 
cess. Assuming the total power of the two pulses is fixed, 
we define the power mismatch ratio R = 1 — Pl/Ph, 
where Pl , Ph represent the lower and higher laser power 
respectively. We use simulation parameters identical to 
that of Fig. IHa) and the results are shown in Fig. [5] 

In Fig. [5fa), the evolution of e nyX starting from the 
birth of helium electrons is plotted. Except for the bal- 
anced case of R = 1, all other curves rise rapidly first 
and then eventually drop to a low level. This reduction 
can be understood from the electron density plot shown 
in Fig. [5fb). After the collision is over, a large number of 
additional helium electrons are liberated in the wake by 
the overlapping of the higher power laser and the wake 
fields. These electrons lead to the emittance growth be- 
cause they are born with transverse positions far off axis. 
However, the majority of these additional electrons do 
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FIG. 5. The impact of laser power imbalance on the trans- 
verse emittance. (a) the transverse normalized emittance vari- 
ation during the trapping process under different laser power 
mismatch ratios, (b) additional helium electrons generated 
on the higher power side for the 30% power mismatch. 



not get trapped (A\I/ > —1) and fall behind leading to 
the drop in the emittance over time. The final emittance 
can be as low as 10 nm as long as the power mismatch is 
less than 30% and with approximately the same charge. 
We note that there is much flexibility in this idea as one 
could experiment with different gas mixtures and collid- 
ing pulse geometries to optimize the tradeoffs in charge 
and emittance. 

In conclusion, we have proposed a new injection scheme 
for PWFA by utilizing transverse colliding laser pulses. 
This scheme can generate extremely low emittance elec- 
tron beams with sufficiently high current. Using simula- 
tions, we demonstrate the possibility of generating beams 
with I p = 0.44 kA (4.6 pC of charge), e njX = 8.5 nm, 
e n,y — 6 nm and longitudinal emittance of 0.06 keV-ps 
leading to a beam brightness greater than 1.7 x 10 19 
A-rad -2 -m -2 , which is three orders of magnitude higher 
than that can be obtained at the LCLS, therefore this 
scheme may have important impact on future FEL re- 
search. 
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